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A computationally efficient, self-consistent complex scaling approach to calculating characteristics 
of excitons in an external electric field in quantum wells is introduced. The method allows one 
to extract the exciton resonance position as well as the field-induced broadening for the exciton 
resonance. For the case of strong confinement the trial function is represented in factorized form. 
The corresponding coupled self-consistent equations, which include the effective complex potentials, 
are obtained. The method is applied to the shallow quantum well. It is shown that in this case the 
real part of the effective exciton potential is insensitive to changes of external electric field up to 
the ionization threshold, while the imaginary part has non-analytical field dependence and small for 
moderate electric fields. This allows one to express the exciton quasi-energy at some field through 
the renormalized expression for the zero-field bound state. 

PACS numbers: 71.35.-y, 73.21.Fg, 78.67.De. 

I. INTRODUCTION 

The quantum-confined Stark effect (QCSE) has attracted a lot of attention since its discover}*!^ due to its diverse 
optoelectronics applications (see for example Ref.3). In quantum well (QW) structures with a layer thickness less 
than the effective exciton Bohr radius, an electron and a hole are forced to orbit closer to each other, and the binding 
energy of the electron-hole pair increases by a factor between two and three. This additional stability makes exciton 
resonances in QW's observable even at room temperatures. More important, these resonances are still observable in 
a uniform electric field perpendicular to the growth direction, which is as large as 50 times the classical ionization 
field for a three-dimensional exciton. The latter also produces a relatively large red shift in the absorption peak up 
to 2.5 times the binding energy of the exciton, while the width of the exciton resonance grows only moderately. Such 
a behavior is opposite to the properties of 3-d excitons or QW excitons in a field parallel to the plane of the QW. 1 ' 2 
The origin of this phenomenon stems from the fact that the electric field induced decay time for QW excitons is 
determined not by its binding energy (~ 10 meV in GaAs/AlGaAs quantum wells), but rather by the quantum well 
barriers' heights (~ 300 - 400 meV). 

Despite the fact that the QCSE was discovered almost 20 years ago and qualitatively is well understood, there 
is no complete quantitative solution for the problem so far. Usually, it is treated numerically with different degrees 
of approximation. The conventional approach consists in the separation of the problem into two consecutive steps. 
The first step in modelling QCSE is to omit the Coulomb interaction of an electron-hole pair and to find single 
particle resonance positions for an electron and a hole separately. In moderate fields (up to 10 5 V/cm) and non- 
shallow quantum wells (quantum wells with more than one discrete level inside) the lowest-lying resonances are so 
sharp that they are often approximately treated as bound states, neglecting the possibility of tunnelling. These simple 
methods include either an infinite well approximation* 2 * 2 **^* where the effect of penetration of the electron (hole) wave 
functions under the barrier is taken into account by introducing an effective QW width^S** or variational calculations 
for finite potential wells 4,7 , which produce only the real part of the exciton quasi-energy. Within a weak-field regime, 
the Stark shift is a small second-order effect in the perturbation theory sense. The absolute value of a shift decreases 
with a decrease of the QW width. Achieving large shifts requires strong fields, where the perturbation theory fails. In 
this regime a significant broadening of the states appears due to the field induced tunnelling. Therefore, it is necessary 
to devise non-perturbative methods for the correct determination of the exciton resonance energies and tunnelling 
life-times in these systems in the presence of a strong field. 

There are several ways of approaching this problem mathematically. The single-particle quasi-bound states in 
QW's in the presence of a uniform electric field have been studied with the help of the time-independent Schrodinger 
equation using Airy functions with real and complex argumentsa^iS, stabilization technique*****^***, finite difference 
method*** and complex coordinate method**^. Another group of methods utilizes a dynamical approach**^******- where 
the time evolution of a tunnelling wave packet is simulated by numerical integration of the time-dependent Schrodinger 
equation. 

It is worth pointing out that all of the above mentioned calculations of decaying quasi-bound states are made in 
the single-particle approximation without consideration of the Coulomb attraction between the electron and the hole. 
Excitonic effects are incorporated at the next step of the calculations, which usually take into account only a real part 
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of the exciton binding energy shift, and are performed with the help of variational methods2^Li^. 

The goal of this paper is twofold. First, we develop an alternative method of studying QCSE that allows us to 
calculate both the position of the exciton resonances and their field-induced tunnelling life-times on equal footing. 
This method is a generalization of a self-consistent approach developed recently by the authors^ for the exciton in QW 
without the electric field. Second, we introduce the model of QCSE in shallow quantum wells (SQW), which allows 
us to obtain analytical results for the resonance positions and tunnelling widths as well as for the modified exciton 
"ground state" wave function. In this model we approximate SQW by a delta- functional potential. Calculations are 
significantly simplified in this model due to the fact that the single-particle one-dimensional problem for a particle in 
a delta-potential and in an applied uniform electric field can be solved exactly. 

The paper is organized as follows. In Sec. [HI we generalize the self-consistent approach of Refll9l for the presence 
of the electric field, and derive main self-consistent field equations for complex resonance energy positions for QCSE, 
using the method of complex scaling. In Sec. II I II we apply our approach to the model of SQW for calculations of 
the effective exciton potential. One of the results of these calculations is the phase diagram representing regions of 
exciton's stability in strong electric fields for different well's widths and/or barrier potential. Concluding remarks are 
given in Sec. IIVI 



II. COMPLEX SCALING SELF-CONSISTENT APPROACH FOR QCSE 



One of the goals of this paper is to develop a new method for calculating excitonic characteristics in QW in the 
presence of the electric field due to QCSE. In order to present the method clearly we consider the simplest model 
of QW excitons. In particular, we assume that both conduction and valence bands are non-degenerate, and that 
they both have an isotropic parabolic dispersion characterized by the masses m e and rrih, respectively. Thus, we 
neglect valence-band mixingSS, anisotropy of the system, and dielectriaSI and effective mass mismatch between the 
well and barrier materials. All these effects, which are very important for realistic computations can be, however, 
easily incorporated into our approach. Throughout the paper we use effective atomic units (a.u.), which means 
that all distances are measured in the units of the effective Bohr radius clb — h 2 e/^,*e 2 , energies in the units of 
/i*e 4 /?i 2 e 2 = 2 Ry, and masses in the units of reduced electron-hole mass fi* , where 1/ fi* = 1/m* + \jm* h , and e is an 
average background dielectric constant of the well and barrier materials. In these notations TO e ,ft — m* h /fJ-*, where 
m* h are effective masses of an electron and a heavy hole. 

Within the approximations explained above, and after the standard procedure of excluding the center-of-mass of 
the perpendicular motion in the plane of the layers^, the excitonic Hamiltonian is given by 



H = E g +H e + H h +K r + V reh , 

H e ( Ze ) : 



(1) 
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where E g is a gap energy, z is the growth direction, r measures a relative electron-hole distance in the transverse 
direction r = \J (x e — Xh) 2 + {y e ~ Vh) 2 ■, are the quantum well confining potential in z direction for the electron 
and the hole, respectively, and F — \e\E measures the strength of the electric field. We also assume that the quasi- 
bound exciton state is independent of the angle in the xy plane, and exclude the corresponding term from the kinetic 
energy of the relative motion K r . 

Strictly speaking, because of the possibility of tunnelling, there are no bound states in a system with uniform electric 
field, and the energy spectrum of Hamiltonian Q is continuous: E G (—00,00). This means that the respective 
wave functions do not vanish at infinity and, therefore, are not normalizable in L 2 sense. Then a proper approach 
to the problem is to consider a scattering problem. The scattering matrix S contains all the physical information 
about the system. For example, resonances correspond to the poles of the 5-matrix in the lower half of the complex 
energy plane. It is well known that these complex energies can also be found by solving an eigenvalue problem for 
the time-independent Schrodinger equation, when the wave functions are chosen to have only outgoing components, 
but do not have an incident wave (so called Siegert or Gamow wave functions). Formally, after matching the usual 
"bound-state" boundary conditions one can obtain the correct complex eigenvalues. However, the eigenfunctions 
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corresponding to these complex eigenvalues exponentially diverge. This circumstance prohibits an application of a 
standard variational principle to calculations of approximate values of the quasi-bound states for QCSE. 

The complex coordinate (complex scaling) approach (see for review Refs. El2l is devised to circumvent this 
problem. The idea consists in applying the coordinate transformation R — * Rexp(^7) to all coordinates in the 
Hamiltonian ft}. Such a transformation makes the new Hamiltonian non-Hermitian [H — ► H(-y)] and changes its 
spectrum respectively. The continuous spectrum moves^ to — ioo. In addition, a discrete set of complex points 
lying below the real energy appears. These complex quasi-energies correspond to the tunnelling resonances under 
consideration. It is important that the eigcnfunctions ^ 7 (R) corresponding to these complex eigenvalues W = E—iT/2 
are square-integrable, and, therefore, can be calculated by means of the variational principle, which requires the 
stationarity of the functional 

W [*] = J * 7 -H"( 7 )* 7 dV I J * 7 * 7 dV. (2) 

It should be noted, however, that the biorthogonality normalization (without complex conjugation of the left eigenvec- 
tor) is used in the definition of the energy functional instead of the standard scalar product. The complex quasi-energies 
should also be stationary with respect to changes in rotation angle 7. The method of complex scaling and methods 
similar to it are widely and successfully used in studies of atomic and molecular auto-decaying states? 2,23 . 

The new variational principle, Eq. (J2J , for the rotated non-Hermitian Hamiltonian can be used for extension of the 
self-consistent method, recently developed by the authors^, to the problem of QCSE. According to this method, we 
construct an approximate exciton wave function ^E' 7 (z e , Zh, r) with the help of the unknown functions ipi,ip2, ■ . ., where 
each function ipk depends on a fewer number of variables than the total wave function. Considering variations of these 
functions independently, from the variational principle, Eq. J2J, we obtain coupled integro-differential equations for 
ipk- The resulting coupled equations can be solved by successive approximations. The convergence of this procedure 
allows for estimating to what degree a given functional dependence of the trial function is close to the exact wave 
function. In the case of zero electric field and sufficiently strong localization (the exciton "z-size" is smaller than 
its three-dimensional Bohr radius), we founcU^ that representation of the trial function in the form of a product of 
functions of different coordinates gives an excellent approximation for exciton energies even after only one iteration. 
We expect, therefore, that a similar form will give accurate results for quasi-bound states in the presence of the electric 
field, at least, when it is not very large. Thus, we choose the trial function for the Hamiltonian in the form of the 
product: 

^trial(r,Z e ,Zh) = -tjj(r)Xe(Ze)Xh(zh)- (3) 

Assuming normalization of every function in this product, we substitute function Eq. in Eq. I|2|l. vary each function 
in a product separately, and obtain the system of coupled integro-differential equations 

[e- 2 ^K r {r) + e-^V r (r)]i;(r) = Wxi>(r), (4) 
[H e (e^z e ) + e'^V e (z e )] X e(z e ) = W eXe (z e ), (5) 
[H h (e^z h ) + e^V h (z h )] X h(z h ) = W hXh (z h ), (6) 

where the following notations for effective potentials are introduced: 

V r (r) = (XeXh\~l/R\XeXh), (7) 
V e , h {z e , h ) = (^Xh,e\-l/R\^Xh,e)- (8) 

The angle brackets of biorthogonal inner product imply that the integration of the Coulomb potential with corre- 
sponding wave functions is carried out over two of three independent variables. 

Solving the system of equations @ (JSJ we obtain the best approximation for the entire wave function in the form 
of a product The corresponding value of the total quasi-bound energy W — E — iT /2 is given by Eq. J3J), which 
can be rewritten in the form 

W = (* = W e + W h + W x - (Xe\e-^V e \ Xe ) - (Xh\e- n V h \ Xh ). (9) 

Eigenfunctions and corresponding complex eigenvalues of Eqs. l|4" )l — can be found with the help of successive 
iterations. As in the case of zero electric field, the zero order solution is obtained after neglecting the effective Coulomb 

terms [^^=0] in single-particle Eqs. JSJ and ©: 

H e , h (e^z e , h ) X f h {z) = W^x^(z). (10) 
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With the help of functions x e L a zer0 or der approximation for the effective potential V r (r) is calculated and sub- 
stituted into Eq. Q . This effective potential has the meaning of a two-dimensional effective electron- hole interaction 
induced by collective effects of an attractive Coulomb interaction, quantum well potentials in the z direction, and an 
applied uniform electric field. The first quasi-bound state for this effective two-dimensional exciton is found from the 
following equation 

^(°)(r) =w£ t V° ) (r). (11) 

The eigenfunction (r) is then substituted into Eqs. JSJ to calculate a new approximation V^ h (z e h) for the effective 
potentials. This process can be continued until the potentials are self-consistent to a high order of accuracy, i.e. until 
the condition 

{i>\VM) « (Xe\V e \ X e) « (Xh\V h \ X h) (12) 

is fulfilled. 

The main advantage of the approach described above is that it allows one to calculate, in principle, not only 
field-induced single-particle widths T e ^ but the exciton width Tx as well, which describes a renormalization of the 
electron-hole pair life-time by the effective interaction. For moderate fields this width is a non-analytical function 
of F, and is small compared to T e ^- The exact calculation of this contribution is possible only numerically. The 
details of such calculations for particular semiconductor structures will be presented elsewhere. Below we consider 
the properties of the solution of Eqs. Q— |JBJ obtained after the first iteration of the self-consistent method for the 
case of a shallow quantum well, where some analytical results can be derived. 

It should be stressed that the most important feature of the complex scaling method is the biorthogonality normal- 
ization. The other feature, the coordinate rotation itself R — > Rcxp(i7), is necessary to guarantee the convergence 
of the integrals 1)718(1 . Although, this rotation is important to negate the asymptotic divergence of the exact eigen- 
functions of the quasi-bound state, it was found^ that in the case of sharp resonances (when the width T is small in 
comparison with the real part of the energy) the optimal angle of rotation 7 is close to zero. Moreover, the final answer 
for quasi-energy is not very sensitive to the change of 7. Therefore, for approximate calculations we will fix 7 — 
and choose the trial functions to be square integrable from the very beginning. However, we will permit variational 
parameters to be complex. The latter assumption allows us to extract the resonance width in the problem, while still 
keeping the problem tractable. We will follow this procedure in the subsequent sections. 



e- 2 ^K r + e-^V ( "'(r) 



III. ^-FUNCTIONAL MODEL FOR A SHALLOW QUANTUM WELL IN AN ELECTRIC FIELD 

A zero approximation of a self-consistent field approach for the electron (hole) part of the wave function in the z 
direction, Eqs. (|10|l . is equivalent to the standard problem of a single-particle quasi-bound state in QW. Therefore, 
its solution can be obtained with the help of any method described p r eviously&2i±2iiAii2iiiLi4ii5ii§iiI. Let us first solve 
the problem exactly using Airy functions with complex arguments. 

We define a shallow quantum well as a well in which only one bound state exists for both electrons and holes. It 
is well-known that, if a well's width, L, tends to zero, a shallow quantum well can be successfully described by a 
5-functional well 

V(z) = U-aS(z), (13) 

where a = UL is characterizes the <5-potential strength and U is the barrier height. For convenience, the first constant 
term U on the left hand side of Eq. I|13|) will be omitted hereafter. In terms of the total Hamiltonian, Eq. (JTJ, it 
means that the energy band gap, E g , is the barrier's energy band gap: Eg 3,1 ' = E™ 011 + U e + Uh- The delta- functional 
approximation is directly applicable only either to very narrow QWi2*2^ or to wells with a small band-gap offset^ 
i.e. when the well width and/or the band offsets are very small so that the carrier wave functions are mostly in the 
barrier region. However, the range of its applicability can be extended^ up to the moment of the appearance of the 
second bound state by introducing the effective strength of the potential 



2U .... 
a = UL cS =\ — ? ( 14 ) 



where x\ is the solution of the equation 



x = 1 — 2/7rarcsin = . (15) 

\V2mt/L 2 / 
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L e g is chosen in such a way that the energy of the bound state in the 5-potential matches the ground state energy 
of the finite well problem^. In some sense, the model of the ^-functional QW is complementary to the model of 
an effective infinite quantum well 2 (EIQW), which is used to approximate finite QW with large widths (or barrier 
heights) when the number of levels in a well is large. Indeed, the more discrete levels exist in the QW the better 
the EIQW model works for the ground state, but it fails when the well has only one discrete level. On the other 
hand, the ^-functional QW is not applicable for quantum wells with more than one level. For typical parameters in 
Alo.3Gao.7As/GaAs structures the range of applicability of 5-functional SQW extends up to a well width L w 40A 

The advantages of the ^-functional SQW model become especially obvious in the case when a uniform electric field 
is applied. The simplicity of the model allows one to derive analytically both the electric field induced energy shift, 
and the tunnelling life-time. Let us first find the complex energy and the corresponding wave function of the single 
one-particle quasi-stationary state for the electron (the case of a hole is obtained by substituting — F for F). The 
respective one-dimensional Schrodinger equation is 



1 d 2 

ao(z) — Fz 



2m dz 2 

It is convenient to introduce the dimensionless variables for the field and the wave vector, 



X (z)=E X (z). (16) 



/ = 2mF/«g, k = y/-2mE/K , (17) 

where kq — ma is the ground state wave vector for an unperturbed ^-potential problem. After performing the 
standard coordinate transformation 



y = K,r" (1 



—) 

Kof J 



(18) 



we obtain 



+ vx(v) = -2.r 1/3 <% + y )x(y), (19) 

where j/o = k 2 /~ 2 / 3 . Everywhere except the point z = (or equivalently, y = — yo), Eq. I|19f) coincides with the Airy 
equation^ In order to find the quasi-bound state we have to use boundary conditions at infinity, which require that 
the solution does not contain incoming waves. It means that we construct our solution in the form of an evanescent 
wave to the left, and propagating outgoing wave to the right: 

Ci(K)[Ai(-y)-iBi(-y)], y -> +00 -^L__ exp (z2/3y 3 / 2 ) z>0 



C 2 (K)M(-y), ^ 2 exp (-2/3|y| 3 / 2 ) 



(20) 
z < 



Note that propagating asymptotic for the right part of the wave function is valid for z 3> k 2 /kq/, i.e. far away from 
the origin z = 0. Near the origin, the wave function decreases exponentially for both z < and z > 0. For small fields 
the evanescent asymptotic can be used to describe the behavior of the wave function at both sides of the point z = 0. 
The unknown constants C\ <x and complex energy are found after matching the boundary conditions at y = —yo- For 
the normalized energy k 2 we obtain the following equation: 



It has one complex solution 



« 2 = £*(/)+M/)- (22) 



The real and the imaginary parts numerically calculated from Eq. 1)21(1 are shown in Figs. Q El The notion of the 
quasi-bound state has sense, when the width ej is small in comparison with the real part of the energy. This happens 
at small fields, i.e. when the condition J 1 / 3 < 1 is fulfilled. In this case the right hand side of the Eq. (|21ll is small, 
and the real part of the argument of the Airy functions is large: Re(K 2 / -2 / 3 ) ^> 1. With the help of asymptotic 
expansions of the Airy functions, 26 Eq. I|21|) can be reduced to 

5 f 2 1155 f i ( 4k\ / 5/ \ 

re=H —x H 4^ H h-cxp 1 — I . 23) 

32 k 2048 k 2 V J \ 24k J V ' 
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FIG. 1: The real part er of the normalized energy re 2 as a function of normalized electric field / = 2mF/t^ ) . The solid line 
represents the exact solution of Eq. 12 U . The dashed line shows the quadratic behavior of an approximate solution, Eq. 12 11 
for small /. The dot-dashed line is the real part for the energy, Eq. 12911 . obtained with the help of a variational method with 
a trial function 12711 . See below in text. 



Solving this equation by iterations, we obtain after the first successive approximation 



1 + ^/ 2 , eiwexpf-^V (24) 



16 J ' V 3// 

For comparison we present these results in Figs. |21 along with the results of the exact numerical solution of Eq.l|2ip. 
The first two terms in the expansion of the wave function l|20[l near the origin with respect to the field parameter / 
are 



X(z) « 



/ 

1 + -sign(z)|K z|(l + |re z| 



(25) 



The result for the real part of the energy presented in Eq. (|24|) can also be obtained from the second order of the 
standard quantum mechanical perturbation theory: 



E M=-2mJ o — — - § + fc2 dk, (26) 

where xfldA 2 -) ~ sin(/cz)/y / 7r are the odd wave functions of the unperturbed Hamiltonian. The non-analytical 
exponential dependence of the field induced broadening in Ea.l|24|) is typical for the Stark effect in a hydrogen- like 
atom r oc exp [— (4/3)(Eq/F£)], where Eq and I are the absolute value of energy and the average radius of localization 
for the ground state without the field. In the case of QW excitons, the main contribution to Eq comes from the single 
particle confinement energy, which is much larger than the exciton ground state energy Eq/Ex ~ U 2 L 2 /Ry 2 a 2 B ^> 1. 
Because of this circumstance, QCSE in QW is observable for very large fields. 

In a strong enough electric field, the broadening quickly grows and becomes comparable to the real part of energy. 
In this case the resonance associated with the quasi-bound state disappears, becoming invisible on the background of 
the continuous spectrum. The condition of the existence of the resonance can be approximately formulated as / < 1. 
Since / = 2mF/(mUL e ff) 3 , it is easy to determine the region of parameters where quasi-bound states (and QCSE 
respectively) are observable. In Fig. [21 we present the phase diagram for the region of the existence of the quasi-bound 
states as a function of well width, L, and electric field for a GaAs-Gao.7Alo.3As SQW. As expected, the larger fields 
and smaller widths destroy the localization of the exciton in the quantum well. 

A knowledge of the complex energy and the corresponding wave function, Eq. (|2(J|I . allows us, in principle, to find 
the corrections to the quasi-bound energy due to the effective electron-hole interaction. To this end, one would have 
to perform a coordinate rotation z e> h — > z e .; l exp(i7) to provide the convergence at infinity for the electron and the 



\(x ( Z dd {z)\-Fz\x[°\z))t 
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Normalized field, f=2mF//% 



FIG. 2: The imaginary part ei of the normalized energy re 2 as a function of normalized electric field / = 2mF/K,Q. The solid 
line represents the exact solution of Eq. 12 U . The dashed line shows the exponential non-analytical growth of an approximate 
solution, Eq. (12 H . The dot-dashed line is the imaginary part for the energy, Eq. 1291 . obtained with the help of a variational 
method with a trial function 1271 . The insert shows the behavior of the real and the imaginary parts of the variational parameter 
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for the trial function l)27|l . 




10 20 30 40 

L, (Angstroms) 

FIG. 3: The phase diagram (L, F) for the region, where the quasi-bound states for an electron (upper curve) and a hole 
(lower curve) exist. The data are presented for GaAs-Gao.7Alo.3As SQW. The calculations are based on the following physical 
constants: ml = 0.067mo, m* h = 0.45rao, U e = 340 meV, Uh = 70 meV, e — 13.8. The upper border for widths, L = 40 A, is a 
limit of applicability of 5-functional model. 

hole eigenfunctions (121)1) . to substitute these wave functions into the integral in Eq. (jJJ), and calculate the effective 
complex potential for the electron-hole interaction. 

Unfortunately, the exact forms for the wave functions H20[l lead to fairly cumbersome numerical calculations of 
the integral (JJJ). However, the simple and rather accurate analytical estimates can be obtained with the help of the 
variational wave function of the following form 



Xtnai(z; (3) = \/kq(1 - /3 2 )exp [«o(-|s| 



(27) 
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A function of similar form has been successfully used in the past to model the real part correction to the single particle 
energy in quantum wells with finite width and in the presence of the electric fieldt In our approach, however, we do 
not restrict the variational parameter (3 to be a real number, extending the range of its values on the whole complex 
plane. The only restriction is Re(/3) < 1, which is dictated by the square integrability of the biorthogonal scalar 
product. Applying the variational principle to the Eq. I|16f) we obtain the following equations for the energy 



K *(f3) = l-f + -±r_. (28) 

The corresponding equation for the variational parameter (3 is 

-f/2 + p~f/2(3 2 - 2(3 3 +(3 5 = 0. (29) 

For small fields, the last three terms on the left hand side of Eq. (|29|) can be omitted, and we obtain (3 « //2. The 
corresponding value for energy shift Ejj-lsi / 2 /4, which differs from Eq. (|24|l only by a factor of 1.25. For moderate 
values of the electric field this variational energy is even closer to the exact value than the perturbative result, Eq. 
H24JI . It is more important, however, that Eq. H29[l also allows one to analyze the field induced decay of the electron. 
In order to understand how the electron acquires a finite life-time, one needs to include the third quadratic term in 
Eq. (|29l) . The solution of the modified equation takes the form 



p = r 1 [i - vw^J , (30) 

from which it is seen that when the electric field exceeds a critical value / = f cr = 1, a square-root singularity appears, 
and the corresponding energy acquires an imaginary part. The exact solution of Eq. I|29|l changes the type of singular 
behavior of function 0(f) in the vicinity of / cr , and lowers the critical field value 

U = 0.487. (31) 

The dependences of the real and imaginary parts of variational energy, Eq H28|) . on the normalized electric field, 
obtained from the solution of Eq. (|29() are shown in Figs.2] ari dGl We can see that the variational results for the real 
part of the energy are in very good agreement (especially in the case of small and moderate fields) with the results of 
the exact solution of Eq. I|21|l . The variational method with a simple trial function l|27|) also captures the effect of the 
field-induced decay of excitons characterized by a non-analytical dependence of the decay rate upon the electric field. 
Thus, the inclusion of the complex solutions of Eq. I|29|) for the variational parameter gives substantial physical insight 
into the problem of calculation of the effective exciton potential. We have to stress here that in the above procedure 
we performed only partial optimization of the functional J3J), since we fixed 7 = 0. The variational optimization of Eq. 
1 01) with respect to the rotation angle would give a non-zero optimal angle 7 op t and even better agreement with the 
exact solution of Eq. I|21|l . As we already mentioned, in the case of sharp resonances^ the optimal angle of rotation 
7 is small. Moreover, the final answer for quasi-energy is not very sensitive to the change of 7. Therefore, the full 
optimization gives minimal corrections to the obtained result for the complex energy and variational parameter (3, 
making, however, the further calculations of the exciton effective potential very cumbersome. 

IV. EFFECTIVE POTENTIAL AND BINDING ENERGY CORRECTIONS FOR EXCITON 

The effective potential for a quasi 2D exciton is given by the integral 

V eS (r;^,^,F) = - r dZe r teM^M**)? (32 ) 
J— 00 J -oo \Jr + [z e — Zfi) 

For electron and hole wave functions we will use functions in a form Eq. (|27|1 with corresponding values of the complex 
variational parameter (3. For further calculations it is useful to introduce the following notations: we will denote wave 
numbers, variational parameters, and wave functions, related to the electron (hole) by index 1(2), for instance wave 
numbers = n^' h \ where i = 1 or 2. We also introduce renormalized parameters KiL,R according to 

kil.r = K\ (1 ± Pi) (33) 
K2L.R = t 2 (lTA) (34) 
A = ApmiF/zs?) = fa + 0u, Pa, Pi > 0. (35) 

(36) 
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With these notations the wave functions (|27|l and effective potential, Eq. J7J, can be rewritten as 



KmKiL J exp(K li z), z<0 



V r (r;ni,K 2 ,F) = - 



dzi 



dz 2 



xIlxIl + x\rxIr . xIlxI r + x\ r xIl 



^Jr 2 + (zi - z 2 ) 2 yjr 2 + (zi + z 2 ) 2 
Vi(K 1L ,K 2L ,r) + Vi(KiR,K2R,r) + Vii(KiL,K 2 R,r) + Vh(k 1r , K 2L ,r), 



(38) 



T7(°) 



where the last line introduces notations for each of four terms comprising the integral representation for V r . As we 
had shown before^ the potentials Vjji can be represented with the help of the function 



T a = T(n a r) 



exp(-2K a t)dt n 

/ l , ^ = o [ H o(2K Q r) - Yo(2K a r)] 

Vr + t 2 £ 



(39) 



where Ho is the zero-order Struve function and Yo is the zero-order Neumann or Bessel function of the second kind. 26 
Then 



V/j/(k q ,k^) 



K\lK\rK 2 lK 2 r 

2k\k 2 



T(npr) ± T(n a r) 



Ka ± Kp 



(40) 



and the final result for the effective exciton potential induced by QW confinement and the presence of the electric 
field takes the form 



T7(0) ( K lK2 (l ~ f3 2 )(l ~ f3 2 ) 
V r (r; Kt,K 2 , F) = 



T\l + T 2 l , T\r + T 2 r Til — T 2 r Tir — T 2 l 



UlL + H 2L KlR + K 2 R n 2R - KlL k 2 l - UlR 



Without the electric field Eg. 1411) is reduced to the result obtained in Rcf. 
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V^\r;Ki,K2,0) = - ' 2 ' 2 lK '' 2 2 [K 2 T(Kir) - K\T{n 2 r)] . 

K 2 K-y 

where it was shown that the properties of this potential are governed by a single parameter 



1/1 



1 



d = o — + — ) = 0a(*i)X2 0*2)10*1 - z 2 ) |xi(zi)xi(z 2 )). 



(41) 



(42) 



(43) 



This parameter has a simple physical meaning of the average square of the distance between the electron and the hole 
in the z direction. It controls the transition between small r and large r asymptotics of the zero-field potential 14211 : 



vf \r; Ki,k 2 , 0) 



iln(r/d), r<d, 
1 r > d 



(44) 



The effective zero-field electron-hole interaction for the exciton in the SQW starts from the true logarithmic Coulomb 
potential of a point charge in two dimensions that smoothly transforms at distances r *~ d to the screening potential 
(|44|l with three-dimensional Coulomb tails. 

The understanding of the role played by the parameter d helps to substantially simplify calculations of the ground 
state energy. For example, we found that the variational energy obtained with a single parameter trial function 



2exp(d/A) — . 

V trial = /,/, ^ exp-V^+c^M 
^A(A + 2d) 



(45) 



coincides with more sophisticated numerical approaches with excellent accuracy. The respective expression for the 
ground state energy can be written as a sum of two terms: Ex — Ey (A) + E^\ where 



4 0) (A) = -^ 



1 



1 



Al + 2d/A A 2 



1 - 



(2d/X) 2 E 1 (l,2d/\) exp(2d/A) 
1 + 2d/\ 



(46) 
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FIG. 4: The profile of the real part of the effective self-consistent exciton potential V r (r; ki, K2,F) for different values of the 
electric field (dot-dashed lines). The thick solid line (curve 5) below the other curves represents — 1/r behavior, and is given 
for comparison. The insert is a magnified view of the circled area. The labels for the curves correspond to, curve 1, F — 
V/cm, curve 2, F = 0.4 x 10 5 V/cm, curve 3, F = 0.8 x 10 s V/cm, curve 4, F = 4 x 10 s V/cm. All data are for d/a B = 0.16 
that corresponds L = 20 A finite quantum well in AlGaAs/GaAs materials. 



with Ei (a;) being the exponential integral^, and 

(47) 



2d 
A 2 



1 2d (2d 
- - 0.557— + 0.563 — 

2 A V A 



The first term E^ 1 should be taken for the optimal value of A that brings Eq. 146|) to the minimum. The second 
term, Eq. (|47|l . is a small correction, which is calculated for a given optimal A. 

In order to understand how these results are modified by the electric field we study the properties of the effective 
potential Eq. (|4ip. Figures 01 and |S] show the profiles of the real and the imaginary parts of this potential calculated for 
several different electric fields. For concreteness we chose for calculations parameters typical for GaAs-Gao.7Alo.3As 
SQW with width L = 20A In our calculations we used the following physical constants: m* = 0.067too, mX — 0.45mo, 
U e = 340 meV, Uh = 70 meV, e = 13.8. For these constants we have K\ — 5.68, K2 — 7.31, and d = 0.16. The latter 
value tells us that in a zero field the exciton is strongly localized in the QW. The critical fields F^ r 2 derived from 
Eqs. (|31I17|I are F^ r = 0.82 x 10 5 V/cm and F e cr = 2.57 x 10 5 V/cm. According to our phase diagram, Fig. a field 
that larger than 1.7 x 10 5 V/cm should be considered as strong, since in this case the hole has already a substantial 
probability to tunnel through the barrier, escaping an interaction with the electron. Results presented in Figure 0] 
show that effects of the electric field on the real part of the effective exciton potential are relatively small. It indicates 
that the corrections to the real part of quasi-bound exciton energy due to the field must also be small. The result is 
consistant with the previous calculations^ for larger QW, where a trend to a substantial decrease of the Stark shift 
was observed with narrowing of quantum well widths up to L = 50A 

With increasing applied field the potential changes in the following way: at small fields (see curves 2,3 on the insert) 
both /3s are real and there is no imaginary part at all. As the field is applied, the electron and the hole are pulled 
apart by the field. It weakens the effective attractions and the curve moves upwards, manifesting a Stark shift for the 
binding energy corrections. Curve 3 represents the maximal deflection from the zero-field potential. It corresponds 
to the moment of the first critical field, when the single particle energy for the hole in QW acquires the width and 
its tunnelling becomes possible. At this moment the optical (imaginary) part of the effective potential also appears. 
After this the real part of the potential begins swinging back. The latter can be understood if we have a look at the 
real part of f3 in the insert of Fig- El After the "critical field" is passed, (3r also changes its behavior and even slightly 
decreases with field growth. The latter explains the "effective" increasing of the attraction again. However, since now 
the connection with the continuum is established, the physical interpretation is not so obvious. The imaginary part 
of the effective potential continues to grow, which is clearly seen in Fig. [SJ In order to understand the character of 
this growth, we fixed the radius of the potential and investigated the field dependence of the potential strength at this 
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FIG. 5: The profile of the imaginary part of the effective self-consistent exciton potential V r (r; F) for different values of 
electric field (dot-dashed lines). The upper solid line (curve 4) represents 1/r behavior, and is given for comparison. The labels 
for the curves correspond to, curve 1, F = 10 s V/cm, curve 2, F = 2. x 10 5 V/cm, curve 3, F = 4 x 10 5 V/cm. All data are 
for djas = 0.16 that corresponds to L = 20A finite quantum well in AlGaAs/GaAs materials. 




12 3 4 
F (10 5 V/cm) 



FIG. 6: The dependence on electric field of an imaginary part of the effective potential V r (r; F) at fixed value of radius 
r/a B = 0.1 Data are for L = 20l AlGaAs/GaAs SQW. 



fixed r. Fig. Ogives corresponding results for this dependence. As we can see, it follows the behavior of the imaginary 
parts of /?i 2 (compare with Fig.EJl. The second leap corresponds to the second critical field, when a new scattering 
channel, namely, a possibility of electron tunnelling, is open. 

The smallness of the correction to the real part of the effective potential suggests that we can try to express this 
effective complex potential at finite field through the zero-field potential, renormalizing some of its parameters. To 
do so, let us look at the asymptotic behavior of the potential Eq. 1)41(1. One can show that the first two terms in the 
asymptotic expansion are 

V ( °\r; Kl ,K 2 ,F) a -1/r + J 2 /(2r 3 ). (48) 
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They are the same forms as for zero-field potential [see Eq. 1)4 4[l]. but with 

1 1 + 3/3 2 1 1 + 3/? 2 4/3!/? 2 



(49) 



It is also easy to see that this new parameter is the average square of the distance between the electron and the hole 
in the z direction, when the electric field applied: 



kXv{z{)xi{zi)\{zi - z 2 f\x\{zi)x\{z2)) ■ 



(50) 



This parameter is complex now, and its imaginary part provides the coefficient for the asymptotic of the imaginary 
part of the effective exciton potential. With this finding, it is easy to understand how we should renormalize the wave 
numbers «i a to adjust corrections in binding energy, induced by electric field. The easiest way is the following: 



1 " 
hi 
1 : 

K 2 



1 1+3/? 2 

■2\2 



k\ (1 
1 1 



01) 
Z01 



«1K2(1 ~ 0l){l 
20102 



01) 



(51) 

nl(l-0lY ' n 1 K 2 (l~0l)(l-0iy (52) 

Then, for small applied fields the potential, Eq. (|41|) coincides with the modified potential for zero field Eq. (|42H up 
to the third degree in the electric field expansion, 



vf\r]Ki,K2,F) 



7T(0), _ s 

V r (r;Ki,K 2 ,0). 



(53) 



Therefore, the expressions ()46I47|I for exciton binding energy in zero field with renormalized parameter d 2 are valid 
for non-zero electric fields. 



V. CONCLUSIONS 



We developed a self-consistent complex scaling approach for calculations of field-induced complex energies of single 
particle states as well as of the exciton quasi-energy in QWs. This technique is a generalization of the self-consistent 
approach recently developed by the authors^ for the exciton in a QW without the presence of the electric field. For 
the case of strong confinement, we represented the trial function in factorized form and obtained the corresponding 
coupled self-consistent equations which include the effective potentials (|7I8|) . These effective potentials are 

results of averaging over z coordinates of the Coulomb interaction and the quantum well potential. We would like 
to note, that, although the resulting equations are obtained for a separable wave function, the approach itself has a 
broader range of applicability, and can be employed for different forms of trial functions. 

We further applied our approach to the case of a shallow quantum well. To this goal we developed the model of 
QCSE in a ^-functional quantum well, which allows us to obtain analytical results for the resonance positions and 
widths and the behavior of the modified exciton "ground state" wave function. This occurs because the single-particle 
one-dimensional problem in delta-potential and applied uniform electric field can be solved exactly. We also compared 
the exact results for single-particle complex energies with the complex scaling variational method results for a trial 
function in a simple form of two decaying exponents with different coefficients of decay. The obtained energy has 
both the real and the imaginary parts, and their values are in a very good agreement with the exact solution for 
quasi-energy within a wide range of the parameters. The non-zero imaginary part gives additional physical insight 
into the problem. The simple variational functions for the electron and the hole allow us to calculate the effective 
exciton potential. This potential has real and imaginary (optical) parts. We found that the real part of the potential 
is insensitive to changes of external electric field up to the QCSE ionization threshold. This allows us to express 
exciton quasi-bound state at some field through the renormalized expression for the binding energy in zero field. 
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